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We study all possible observables in B ^ with new physics (NP), including new vector, 

scalar and tensor interactions, and investigate the prospects of extracting NP Wilson coefficients 
with optimal observables. Analysis of the full integrated branching fractions of B ^ 
show that the overall sensitivity of the observables of B ^ DrVr is more towards the scalar current, 
whereas the bin-by-bin analysis of distribution of the differential branching fraction points to 
regions of q^ sensitive to tensor interactions. Interestingly, the observables in B ^ D*Ti>r are more 
sensitive to tensor interactions, and bin-by-bin analysis of this mode shows the distinct regions of 
q^ sensitive to vector or scalar interactions. In addition to that, the r polarisation asymmetry is 
found to be more sensitive to NP compared to the other observables, in both decay modes. 

PACS numbers: 


I. INTRODUCTION 

Measurements of branching fractions and other related observables in semileptonic decays of B meson to r can be 
interesting for an indirect probe of NP. Earlier measurements on = BIZ{B —> ^^TVr)/Bll[B ^ by 

Belle [1] and BABAR collaborations [2] have shown some deviations from their Standard Model(SM) predictions [3], 
indicating a possible signature of NP in 6 ^ crv^ transitions. Several authors have tried to explain the observation in 
various NP scenarios mi, as well as in a model-independent way 019]. In order to distinguish between the possible 
signatures of NP, the study of NP in distribution of differential branching fractions in B ^ various 

correlations among r forward-backward asymmetry and r polarization asymmetry has been examined in the Ref. 
laiiQiin]. Recently, a 2.la deviation in the measurement of R(D*) has been reported by LHCb collaboration [12]; 
Belle collaboration has also announced their most recent results on and measured values are consistent with 

the SM within error bars m- 

So far, the constraints on the new couplings are obtained assuming their presence one at a time 1311 ]. If we consider 
all the interactions together, then it will be an impossible task to extract all the couplings from a single measurement. 
However, if one can reduce the number of coupling parameters by imposing certain constraints on the full set of 
parameters, only then it is possible to obtain meaningful errors on the couplings, although the information lost due to 
various assumptions cannot be retrieved. Therefore, it will be useful to have independent couplings, parametrised in 
such a way that the measured errors on different parameters are uncorrelated. On the other hand, it is not necessary 
for a particular observable to have equal sensitivity to different types of NP operators. Therefore, it is useful to know 
how an observable can be optimised to guarantee the maximal sensitivity to a particular type of NP interaction, 
which in turn will help us select observables suitable for the extraction of a particular type of coupling. Hence, from a 
phenomenological point of view, it is important to find out the significance of different types of NP interaction, to an 
observable. To achieve this goal, we use the optimal-observable analysis using the invariant mass squared g^, of the 
lepton-neutrino system. We construct the optimal observable to identify the NP structure that can be best extracted 
from a particular observable, with reasonable statistics. It also provides a deeper understanding of the sensitivity 
that can be best obtained by any method, for a certain process. This technique has been widely used in collider 
phenomenology [SHIS]. 

In this article, we analyse the g^ distribution of the differential branching fractions, the r polarisation 

asymmetries, forward backward asymmetries of the decay B and D* polarisation asymmetry. We include 

all possible non-standard four-fermi effective interactions of the lowest dimension, and estimate the expected statistical 
uncertainties in the extraction of various NP Wilson coefficients that can contribute to B ^ 
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II. METHODOLOGY 


The most general effective Hamiltonian describing the b crUj- transitions with all possible four-fermi operators 
in the lowest dimension is given by [8] 


^e// = 


(1 + Cvi)Ovi + CV 2 OV 2 


+ CsiOs:^ + CS 2 OS 2 + CtOt 


( 1 ) 


where the operator basis is defined as 

Ovi = (cL7^^L)(rL7M^rL), 0^2 = {cR-f^bR){fL^^Vrh), 

Osx = {cLbR){fRUrL), Os2 = {cRbL){fRUrh), 

Ot = {cRCF^'"bL){fRCF^yUrL), ( 2 ) 


and the corresponding Wilson coefficients are given by Cw{W = Vi,V 2 , Si, S 2 ,T). In this basis, neutrinos are assumed 
to be left handed. Our main focus is on the distribution of differential decay rate dT/dq^ in B ^ The 

complete expressions are given in ref. [5]. 

As mentioned earlier, the optimal-observable analysis is a technique to systematically estimate the statistical un¬ 
certainties of the measurable parameters while extracting them from some observable. Elaborate discussions on this 
technique can be found in references p!5Hl8] . In order to apply this technique to B ^ h is necessary to 

express the distribution of the differential decay rate as 


dr{B TUr) 

dq^ 




(3) 


where Qs are functions of C^s. The theoretical expressions for Qs, along with the /i(g^)s, can be extracted from a 
direct comparison between the similar terms on both sides of eq. The coefficients Q, relevant for the bra nchin g 
fractions in B ^ are given in Tableand the corresponding fi{q‘^)s are given in the Appendix (Table XIV). 


^^Obs 

Ci 

dBjdq^ in H ^ Dti^t 

dBjdq^ in H ^ D*rVr 

Cl 

1 -h Cvi + Cv2 ^ 

1 + Cvi + \Cv2 ^ 

02 

\Csi+Cs,V 

Re[{l + Cvi)C*v,] 

03 

\CtV 

\Csi-Cs,V 

04 

Re[{l + Cvi + Cv 2 ){Cs^ + <^52)] 

\Ct\‘‘ 

05 

Re[[\ + Cvi + Cv^)Cl] 

Re[(l + Cvi-Cv,){C*s, - 05 ,)] 

Ce 

- 

Re[{l + Cvi)C^\ 

C-r 

- 

Re[Cv.,Cl] 


TABLE I: as defined in eq.(|^. The observable P^{q‘^) contains the same set of CiS. 


The goal of this technique is to extract Qs, which can be done by defining suitable weighting functions Wi{q^) such 
as Ci = / Wi{q^){dV/dq^)dq^. In general various choices of re^s are possible. However, there is a unique choice for 
which the resulting error in the extraction of Ci is minimized and these functions are given by 


wM^) = X 

3 


dT/dq^ ’ 


where Xij is the inverse of Mij which is defined as 


Mij — 


/ 


fsM{q^) 


(4) 


(5) 


^ CiS are minimised in a sense that the whole covariance matrix is at a stationary point in terms of varying the functional forms of Wi{q^) 
while maintaining f Wi{q^)fk{q^) = Sik- 
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In the above expression, fsM{<f) can be obtained from eq.(|^ by setting Cw = 0, while Csm = 1- Hence, using eqs. 
Q, and ^ the statistical uncertainties in C* extracted from the branching fractions can be obtained as [THl ITT] 


\SCi 


I XuBjB ^ D(*)T Ur) 

N,in 


exp 


Xii 


crpCeff ’ 


( 6 ) 


where = (l/F) J dq^dT/dq^ is the total branching fraction in the decay B —> with V as the total decay 

width. Nsig is the total number of events. As given in eq.(|^, the errors are also related to the production cross 
section ap ( = and the effective luminosity Ceff = Cint^s^ where Cint and are the 

intcCTated luminosity and reconstruction efficiency respectively ^ . The above-mentioned method, and the equations 
like^^ and ([^, can be generalised for any other observables in B decay. 

Since the data is consistent with the SM, if there is NP in B decays, the effect is expected to be small 

compared to their SM counterpart. The earlier model independent analysis [8], which is based on data by BABAR [2], 
shows that zero value of the new Wilson coefficients are consistent with the data. Therefore, we choose our starting 
point as Cw = 0 and find out errors in the extraction of those coefficients around that point. In addition to that, we 
assume that the error on Ci could be captured sufficiently well by just the leading-order terms. 

We focus on the following observables: 

• The branching fractions, obtained by integrating the differential branching fractions over the full q^ region, 
normalised by the full q^ integrated branching fraction = B{B 


with 




Rd(*) (q^) 


1 dB{B ^ D(*'^TVr) 


(7) 

(8) 


• r polarisation asymmetry which we defined as P^'' \q^) = Pr{q‘^)RDM{Q'^)j where 


Pr{q^) = 


drx=i/2/dq^ - dT^=_^/2/dq^ 

dTx=i/2/dq^ + dTx=-ii2ldq‘^' 


(9) 


• r forward-backward asymmetry (q"^) = AFB{q^)RDi*){d^)^ where 

,2^ _ fo 0 dcosd-f°^ g dcos^ 


■Apsiq^) = 


hi dq-fcose dcosO 


bejq'^) 

dT/dq 


,2 ’ 


( 10 ) 


where 0 is the angle that r makes with the B in the rest frame of ru. The expressions for bo{q^) are given in [5]. 
• T)* longitudinal polarisation asymmetry P^^{q^) = where 


^ _ dq^ V ^ _ 

^{Xn,=o) + ^iXD^=i) + ^{XD^=-iy 


^ As we know that the cross section (Ja^b = craT^j/P, therefore, we can define (Pb^d(*)tu ~ (JpB{B —)■ D^'^^riAr): where ap is the BB 
production cross section. If we redefine our observable as than the errors in Ci can be written as 












he// V (PP^eff 


XiiB{B DiCrurY 

Xsin 


since X'-- = Xu!ap. 
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^^Obs 

Ci 

A§B(g^) 


PD.(g^) 

Cl 

1 + Cvi + Cv2 ^ 

ii + cviV - \cvy 

1 + Cvi — Cv2 ^ 

C 2 

Re[(l + Cvi + Cv 2 ){Cs^ + Cs^)] 

1 + Cvi — Cv2 ^ 

\csi-csy 

Cs 

Re[{l + Cvi + C'vJCJ] 

\CtY 

\Ct\'' 

C 4 

Re[{Csi+Cs2)Cl] 

Re[{\+Cvi-Cv2){C*s, -QJ] 

Re[{l + Cvi - Cv,)(Q, - QJ] 

Cs 

- 

i?e[(l + Cvi)C^] 

Re\{\ + (7vi ~ Cv2)C^\ 

Ce 

- 

Re[Cv2Cy[ 

- 

C 7 

- 

Re[{Csi-Cs2)Cl] 

- 


TABLE II: Expressions of CiS for different observables. 

In the above definitions, the detailed expression for dT/dq^ are taken from ref. [5]. For forward backward asymm etries 
and the D* polarisation, the C^s and the corresponding /i(g^)s are given in the Tables \n\ and in the Appendix XVI 


All the above-mentioned observables are expected to be measured with good statistics in future experiments like 
Belle-II and LHCb. The corresponding errors on Ci can be obtained using the following relation 


where X/, = XuBe. 


\5Ci\ = 


xi 


B( <7p Ceff 


( 12 ) 


III. ANALYSIS 


There are varieties of NP models that can contribute to B and the characteristics of those models 

could be very different. For example, two higgs doublet model (2HDM) has only scalar-type interactions, new gauge 
boson Z' and W' take part only in vector-type interactions, the model with leptoquarks has both the scalar or vector- 
type interactions [ 20 ] the extra dimensional models have tensor interaction in addition to scalar or vector-type 
interactions [ 21 ]. 

In our analysis of the decay B —> h will be hard to estimate the uncertainties in the extractions of Cy^ and 

( 7^2 5 because they cannot be singled out from their SM counterpart (same /^s). The similar argument holds for the 
decay B however, in this decay we can estimate the error in the extraction of ReiCy^)^ see for instance 

Tables |l] and XV where /2 associated with C 2 is different from fi associated with Ci. In Table |ml we list a few 
interesting cases of NP relevant for the observables in B ^ In many cases, we assume Cy = Cy-^ -h Cy^ = 0, 

however, the assumption Cy 7 ^ 0 will lead to the same set of parameters that has to be simultaneously extracted, if it 
is assumed that Cy << 1. Under such conditions Ci can be treated as the Wilson coefficient of the vector operator. 
The different NP cases related to the observables in B are given in Tables IV[ V] and respectively. 

In most cases, we assume Cy^ = 0, though the same set of parameters can be obtained without this assumption if 
Cvi << 1. For r forward-backward asymmetry in B ^ we discuss mostly the cases with Cy^ =0. In such 

cases. Cl = C 2 ^ and therefore we need to merge /i(g^) and / 2 (^^) into f{q^){= fi{q^) + / 2 (^^)) for the analysis. In 
all the other cases, when Cy^ 7 ^ 0, the extracted uncertainties are large. We will discuss only one such interesting 


Cases 

Assumptions 

a 

C, /O, z = 1,..5 

b 

ReiCr) = 0 

c 

Re{Cs) = 0 

d 

Re{Cs) — 0 and ReiCr) — 0 

e 

0 

II 

f 

II 

0 


TABLE III: Cases relevant for RD{q^)^ P^{q‘^) and Here Cs = Csi + Cs 2 ^ and in all the cases Cy^ + Cy^ = 0. 


The numerical values of all the relevant parameters, like the form-factors, various masses and lifetimes are taken 
from ref. [22], and for the analysis we choose the central values of all the form-factors. The errors of the form-factors 


Although the model with scalar-type interactions may also contribute to tensor Wilson coefficients by Fierz reordering. 
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Cases 

Assumptions 

a* 

II 

o 

P 

O 

II 

o 

c* 

o 

II 

(P 

Re{Cs) = 0, Cva = 0 

e* 

Re(Cs) — 0, Ct — 0, Im{Cy 2 ) ~ ^ 

/* 

Cs = 0, Re{CT) = 0, Im{Cv^) = 0 

9* 

Cv 2 = 0, ReiCr) = 0, Re{Cs) = 0 

P 

(7v2 = 0? Ct — 0 

P 

Cy^ =0, Cs = 0 

f 

Re{Cs) = 0, ReiCr) = 0, ImiCv^) = 0 


Cs — 0, Ct — 0 


TABLE IV: Different cases related to RD*{q^) and (q^) in B ^ Here Cs = Csi — Cs2^ and in all the cases 

_ Cy^ — 0. _ 


Cases 

Assumptions 

1 * 

II 

o 

2 * 

Cy2 = 0 , Rc^Ct) — 0 

3* 

Cv 2 = 0, Relcs) = 0 

4* 

Cv 2 = 0, ReiCr) = 0 , Im{Cs) = 0 

5* 

Cv 2 = 0 , Cs = 0 

6 * 

Cv 2 = O 5 Rg-{Cs) — 0 , Rc^Ct) — 0 

7* 

Cv 2 = 0 ? Ct = 0 


TABLE V: Cases relevant in Apsiq^) with Cy-^ — 0. 


Cases 

Assumptions 

A 

Ci^O, i = 1,..,5 

B 

ReiCr) = 0 

C 

Re{Cs) = 0 

D 

Re{Cs) = 0 , ReiCr) = 0 

E 

0 

II 

F 

II 

0 

C 

Re{Cs) = 0 , Ct = 0 

H 

ReiCr) = 0 , Cs = 0 


TABLE VI: NP cases relevant in D* polarisation asymmetry. Here, Cy — Cy-^ — Cy^ — 0. 
are considered while we estimate the additional errors on the extracted coefficients. We choose as benchmark values 
B{B Diu) = 2.32%, B{B DHv) = 5.31%, ap = 1105.63 pb, and Ceff = 1/^“^- 


Decay 

B Dryr 

Cases 

b 

c 

Obs. 

RiD) 

pi? 

RiD) 

pi? 

(5Ci 

0.60 

0.37 

0.60 

0.37 

Jjc^F 

1.03 

0.04 

0.13 

0.08 

6\Ct\'^ 

0.62 

0.70 

1.12 

0.72 

SiReiCs)) 

1.31 

0.06 

- 

- 

SiReiCr)) 

- 

- 

1.15 

0.12 


TABLE VH: Numerical values of the la error on CiS extracted from R(D) and P^. Eor the cases Re{Ci) = 0, 


In Table VIb we list our main results of the un certainties in Ci extracted from the analysis of the R{D) and P^ 


corresponding to different cases listed in Table Ilb while those for and P^ , corresponding to the cases listed 


in Table Ell are given in Table |VIII[ For a given case, we estimate the statistical significance of the simultaneous 
extraction of Qs. The numerical values are given only for parameters relevant to a particular case, while the rest are 
set to zero. 

In some more simplified cases, where the number of non-zero NP parameters are less, we compute the x^, which is 
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(a) Case d {R{D)) 


(b) Case e {R{D)) 


(c) Case / (RiD)) 





0.06 

0 Re(Cr) 


(d) Case d (P^) 

FIG. 1: Surfaces of constant ^ 


(e) Case e (P^) (f) Case / (P^) 

for a few selected cases of different observables in B ^ DtVt 


0.39 


- 0.04 


0.86 


1.14 


ImiCsV 0 


Im(Cr) 0 


- 0.15 


0lm(Cr)2 


- 0.39 


0 Im(Cs) 



(a) Case /* (RiD*)) (b) Case g* iR{D*)) 



Im(C,v) 

0 



Im(Cr) 


Re(Cr) 


(d) Case i* (RiD*)) 


(c) Case /i* iR{D^)) 



Im(Cr) 


Re(Cv2) 


0.97 



Im(Cr)^ 


0 Im(Cs) 


2 



(e) Case/* (P^*) (f) Case g* (Pf*) (g) Case h* (P^*) (h) Case i* (P/**) 

FIG. 2: Surfaces of constant = 1 for a few selected cases of different observables shown in P —?■ D*ti>t- 
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Decay Modes 

B D^Th'r 

Cases 

a* 

b* 

c* 


e* 

Obs. 

ParamsT'''''---^ 

R{D*) 

pR* 

R{D*) 

pR* 

R{D*) 

pR* 

R{D*) 

pR* 

R{D*) 

pR* 

6Ci 

7.22 

13.70 

289.17 

116.82 

28.44 

14.08 

2.01 

0.65 

1.28 

1.25 


- 

- 

629.08 

204.37 

56.25 

29.65 

1.00 

1.73 

3.68 

1.83 

6\Ct\'^ 

4.30 

1.96 

11.86 

4.62 

- 

- 

0.03 

0.04 

- 

- 

SiReiCs)) 

- 

- 

529.3 

191.49 

6.81 

2.20 

- 

- 

- 

- 

5{Re{CT)) 

28.27 

36.92 

36.71 

45.37 

- 

- 

0.35 

0.24 

- 

- 

5{Re{Cv,)) 

47.10 

18.51 

- 

- 

14.21 

7.03 

- 

- 

0.63 

0.63 

6Re[Cv2CT] 

15.70 

17.19 

- 

- 

- 

- 

- 

- 

- 

- 


TABLE VIII: Numerical values of the Icr error on CiS extracted from R{D*) and . For the cases Re{Ci) = 0, 

5\Ci\^ = Silm^iCi)). 

defined as 

X" = 

where, Vij = -—-. (13) 

oe ap Leff 

The C-^s are the seed values, which can be considered as model inputs; as discussed earlier, we choose (7-^ = 0, for 
i 1, and = 1. The = 1 surfaces are perfect ellipsoids in Ci basis, and they indicate the Tier errors in the 
determination of Qs. The constant = 1 surfaces are shown in Figs, [l] and[^ The largest and the smallest values 
in the figures represent Tier errors of corresponding parameters. 


Cases 

b 

c 

d 

e 

f 

|5Ci| 

0.27 

0.27 

0.04 

0.07 

0.07 

|<5Pe(C's)| 

0.40 

- 

- 

0.05 

- 

\5Re{CT)\ 

- 

0.30 

- 

- 

0.04 

\5{Im{Cs)Im{CT))\ 

0.11 

0.11 

0.01 

- 

- 


TABLE IX: Numerical values of \5Ci\ extracted from the r forward-backward asymmetry in B ^ Drvr. 


A. Discussions 


For case a, the uncertainties obtained from 
observable R{D) shows that ^ 


the simultaneous extraction of all the Wilson coefficients from the 


\SC 4 \ \SRe{Cs)\ 

|<5C5| \SRe{CT)\ ^ ’ 5\Ct\^ 


which shows that R{D) is equally sensitive to the real part of Cs and Ct- The above result does not allow a direct 
comparison between the sensitivities to the imaginary part of the coefficients. The results obtained for all the other 
cases are shown in Tab. |vnl and in Figs fTa| [T1^ and [T^ We note that if the Wilson coefficients are purely imaginary 
then R{D) is more sensitive to Im{Cs) compared to Im{CT)- Also, it is important to note that this observable is 
more sensitive to the real part of the coefficients than the imaginary part. Comparing all the different cases considered 
for R{D), it would be difficult to comment on the overall sensitivity of this observable to a particular type of NP 
interaction. However, in the next section we will see that there are distinct regions of which are sensitive to either 
scalar or tensor type interactions. 

On the other hand, the analysis of P^(g^) for case a gives 


|^C4| \6Re{Cs)\ 

IJCsl ^ \5Re{CT)\ ^ ' ’ 6\Ct\^ 


( 15 ) 


4 


The results corresponding to the case a are not shown in the table, because the extracted uncertainties are very large (>> 1). 










































which shows an improvement in sensitivity to Re{Cs) compared to Re{CT)- The results obtained from all the other 


relevant cases are shown in Table [VII| and in Figs. andpT[ which allow a case by case comparison between the 

results obtained from R{D) and Interestingly, the extracted uncertainties are less compared to that extracted in 
R{D). The sensitivity of this observable to tensor interaction is a little less, compared to scalar interaction, but it 
can be extracted wit h unc ertainties less than 1 . 

As shown in Table VIIT and in Fig. the observables like R{D'^) and are more sensitive to \Ct\‘^ compared to 
any other new Wilson coefficients, almost in all the cases \Ct\‘^ and Re{CT) can be extracted with small uncertainties. 
A case by case comparison shows that the above observables are more sensitive to Ct than Cs, and SRe{CT) < 
6 Re{Cv 2 ) t>nt they are of same order. Also, when Ct is purely imaginary, we find 6 Im{CT) ~ 5 Re{Cv 2 )i though 
have little better sensitivity to Im{CT)> Therefore, these observables alone won’t allow us to distinguish the 
contributions from right handed vector current to that of a tensor current. However, in the next section we will see 
that a bin by bin analysis of the distribution of the differential decay rate allows to discriminate the effects of these 
interactions. 

However, when Ct = 0 , both the observables are almost equally sensitive, though 6 Re{Cs) < 5 Re{Cv 2 )i to the 
real part of the vector and scalar Wilson coefficients. In case Cs is purely imaginary, 5 Re{Cv 2 ) SIm{Cs) i-c the 
observables are less sensitive to the imaginary part of Cs compared to the real parts of Cv 2 and Cs- Again, we note 
that the extracted errors on Qs from P^ are smaller than those in 


Cases 

2* 

3* 

4* 

5* 

6* 

7 * 

|5Ci| 

1.30 

2.41 

0.40 

0.27 

0.52 

0.04 

S\Ct\'^ 

0.12 

0.08 

0.04 

0.04 

0.03 

- 

\SRe{Cs)\ 

16.30 

- 

1.72 

- 

- 

0.02 

\SRe{CT)\ 

- 

1.06 

- 

0.12 

- 

- 

\5{ImiCs)ImiCT*))\ 

2.48 

2.41 

- 

- 

0.26 

- 


TABLE X: Numerical values of \SCi\ extracted from r forward backward asymmetry in B ^ D^rUr- 


The results of the analysis of forward-backward asymmetries and I)* polarisation are given in Tables [IX| and 
0 respectively. The forward-backward asymmetry in B —> DtVt is equally sensitive to the scalar and tensor type 
interactions. 

For case 1* in we find 


and 


5Ci 


SC 2 

S\Ct\^ 


24, 


(16) 


5C2 

SReiCr) 


12 , 


(17) 


where Ci and C 2 are the functions of the Wilson coefficients of the vector operators. The approximate forms are 
given by 


Cl ^ 1 + 2i?e(CvJ, C 2 ^ 1 + 2Re{Cv,) - 2i?e(CyJ. (18) 

It indicates that the r forward-backward asymmetry in H is more sensitive to tensor Wilson coefficients 

than to a vector, in particular to Cv 2 • In order to understand it better, we define 

Ci2 = Ci-C2^2i?e(CyJ. (19) 

Therefore, a simple calculation shows that 

5Re{Cv^) ^ 1 5Ci2 _ , . 

5Re{CT) 2 8Re{CT) ^ ^ ’ 

In all the other cases with Cv 2 = O 5 the Afb in B is more sensitive to the tensor interaction compared 

to the scalar. On the other hand the D* polarisation is equally sensitive to the scalar and tensor interactions. 
Therefore, if future data shows large deviations from the SM predictions in all the observables like i?(T)*), and 

1)* polarisation, that can be thought of as an indication of the presence of a new tensor type interaction. On other 
hand, if a deviation is only in R{D'^) and not in the others, that could be an indication of a new vector interaction. 
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Cases 

A 

B 

C 

D 

E 

F 

c 

H 

|(5Ci| 

3.58 

0.74 

1.41 

0.12 

0.05 

1.41 

0.05 

0.1 


16.63 

3.39 

0.55 

0.55 

0.77 

- 

0.54 

- 

(5 Ct 

2.53 

0.09 

0.41 

0.01 

- 

0.41 

- 

0.01 

\me{Cs)\ 

7.05 

1.40 

- 

- 

0.20 

- 

- 

- 

\5Re{CT)\ 

5.07 

- 

1.00 

- 

- 

1.00 

- 

- 


TABLE XI: The results obtained from the analysis of the D* polarisation asymmetry in the decay B D^rUr- 

B. Bin-by-bin analysis 


In general, the sensitivity to various NP interactions may also be dependent. Hence, we analyse the bin-by-bin 
distribution of the differential decay rate of H -> to look for more possibilities, and zoom in to the regions 

of g^, within which the sensitivity to a specific type of new interaction is much larger than most other regions. In 
general the 5CiS extracted from individual bins are very large, therefore in the figs. [^|^we plot 5Ci = ^CilNnorm^ 
where Nnorm is some number used to normalised 5Ci. 


0.4 

0.3 

0.2 

0.1 

0.0 

(a) Case b with C 2 = \Cs\‘^, 
C3 = /m(CT)^ C4 = 

^norm — 


Re{Cs)^ and 
^5 




(b) Case c with C 2 — Im(CsY^ 
Cs = |CTr, C^^ReiCr), and 
Nnorm = 10 ^ 



(c) Case d with C 2 — Im{CsY^ 
Cs = /m(CT)^, and Nnorm = 10^. 



4 5 6 7 8 9 10 11 

qhGeV^/ch 

(d) Case e with C 2 = 

C4 = Re{Cs), and Nnorm = 10^. 

FIG. 3: Selected cases in B 


4 5 6 7 8 9 10 11 

q^(GeV^/C^) 

(e) Case / with C 3 = \Ct\‘^, 

C 5 = Re{CT), and Nnorm = 10^ 
here, Cl — \dCi\l Nnorm • 


The results obtained from the analysis of the g^ distribution of differential decay rate in H ^ DrUr are presented in 
fig-i where the variations of the 5CiS with g^ are shown. The normalised uncertainties in the simultaneous extraction 


of |CtP, Re{CT) and Re{Cs), and their variations with g^ are shown in figs. |3a| an d 3b respectively. On the 

other hand the variations of 


and 6Re{Cs) with g^ when Ct = 0 are shown in figTISm while that for 6\Ct\' 


and 


6Re{CT) when Cs = 0 are shown in fig. ^ We note that in the low g^ region (< 7GeV^/c^) the differential decay 
rate is sensitive to the scalar interaction , the sensitivity to tensor interaction in this region is very weak, whereas 
in the high g^ region (> 7GeV^/c^') it is rather sensitive to the tensor interaction. In case the Wilson coefficients are 
purely imaginary, in all the g^ regions, the decay rate distribution is sensitive more to the scalar interaction (fig. 3c) 
than the others. 

As we noted earlier, R{D*) is equally sensitive to Cv 2 and Ct (case a*), however, the analysis of the differential 
decay rate distributions show that (figs, [^and 4d) it is sensitive to tensor interaction only in the very high and low 


^ In very low regions the distribution is also sensitive to Ci . 
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regions, and it is sensitive to Cv 2 in all the regions except the very low region. In fig. |^the variations of 
5Re{Cv2)i and with the in the case Ct = 0 are shown, we note that in all the g^ regions the 

decay rate is equally sensitive to vector and scalar interactions except in the very low g^ region, where the sensitivity 
to Re{Cs) is better than that to Re{Cv 2 )- We also study the cases when the NP interaction is scalar type. The g^ 
distribution of the extracted errors on the respective parameters is shown in fig. which indicates that the decay 
rate is sensitive to scalar interactions only in the low g^ region. If Cs and Ct are purely imaginary than the g^ 
distribution of the decay rate is sensitive to the tensor interactions in all the g^ regions (fig. 4b). 



(a) Case c* with C2 = Re{Cv2), 
C3 = \Cs\^C5 = Re{Cs),8.nd 
at — in6 



(b) Case f with C2 = Re{Cv2), 
C3 = Im{CT)‘^, C4 — Im{CsY, and 

Nnorm = 10^ 



(c) Case h* with C3 = \Cs\‘^, 

C 5 = Re{Cs), and Nnorm = 10^ 



0.004 



5.5 6.5 7.5 8.5 


(d) Case i* with C 4 — |Ct|^, 
Cq = Re{CT), and Nnorm = 10^ 


(e) Case /c* with C2 — Re{Cv2) and 

Nnorm = 10^ 


FIG. 4 : Selected cases in B ^ D*TUr, here, Ci = \SCi\/Nnorm‘ 

All these studies suggest that we could gain in NP sensitivity if we focus on specific g^ regions, which we may 
lose in the full g^-integrated observables. We note that the sensitivity to a particular type of interaction is limited 
to particular regions of g^. Therefore, the experimental data in the specific regions of g^ could help us a better 
interpretation of the type of NP interactions, which may not be obtainable from g^ integrated observables. As for 
example, if we see large deviations in data only in the very high and very low g^ regions, that can be interpreted as 
due to the presence of a tensor interaction (fig. 4d). On the other hand, if we see deviations only in the low g^ (< 7) 
bins, that could be due to a new scalar interaction (fig. 4c). Finally, if data shows deviations in most of the g^ bin 
except the very low g^, this can be due to the presence of new vector interaction (fig. [i^ . 

In order to explain our point we take the example of the g^ distributions of the measured events in B ^ 
which are shown in Fig.0. The plots are generated using the data given in ref. la EH El]. The predicted g^ 
distributions in SM with the central values of the form factors are shown by black lines, blue dotted lines represent 
the errors in SM. In Fig. (5a), we see that the data is not fully consistent with the SM in the region 8.0 < < 11, 


which is the region where the decay rate distribution is sensitive to tensor interaction, as analysed above in the decay 
B DtUt-. At the moment it is hard to conclude anything, and we have to wait for better statistics. From Fig. (5b) 
we see distinct regions of g^ where the data is not fully consistent with the SM prediction, and our analysis suggests 
that those regions could potentially be very sensitive to NP effects. Again, because of poor statistics it is premature 
to conclude anything. Therefore, the experimental effort should be in gathering more statistics in specific regions of 
g^ potentially sensitive to NP, which may in turn help the clean extraction of NP couplings. 

As mentioned above, the numerical estimates are done with central values of all the relevant parameters taken 
from [22]. Numerical instability of our results could be main source of uncertainty in our estimates. The numerical 
results of which depend on the matrix inversion of are often unstable; even a tiny variation of Mij could 

change Xij significantly. This is why, when we estimate the statistical uncertainties in simultaneous extractions of 
the Wilson coefficients, we allowed only stable solutions. We calculate the selected SCi first to and then to 
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(a) (b) 

FIG. 5: Measured distribution of the events in B ^ DrUr (left) and B D*TUr (right) decays [2]. The black lines 
represent the distribution predicted in the SM, obtained from the respective differential branching fractions. 


(m — 1)^^ decimal places, and obtain and respectively. We consider the results as stable only when 

/SCi^'^^ < 0.01. We checked the stability up to m = 10, and in most of the cases presented above, 
our results are very much stable, and the error due to this is negligible. As we can see from the expressions of SCi, 
the other sources of errors in our estimates are given by the errors in /^(g^), ap , jCeff ? and Bi. It is straight forward 
to estimate the errors due to crp , jCeff , and &£. 


Cases 


SCi 

dCi+ 

SCi- 

±%Err. 


Cl 

0.082 

0.079 

0.087 

5.032 

d 

Im{CsY 

0.013 

0.013 

0.014 

5.032 


Im{CTY 

0.168 

0.160 

0.177 

5.032 


Cl 

0.240 

0.229 

0.253 

5.031 

e 


0.274 

0.261 

0.289 

5.031 


Re{Cs) 

0.354 

0.337 

0.373 

5.031 


TABLE XII: Numerical values of SCiS, and SCi~^ {SCi ) considering +10% (-10%) errors in fi{q^) for few cases of R{D). 

The % error is given by {SCi — SCi^)/SCi. 


However, the estimate due to fi{q‘^)s are not that straight forward since M^^s depend solely upon them. The main 
sources of uncertainties in /^s, including the SM, are the form-factors. Errors due to other parameters, like CKM 


element etc, are canceled in the ratios. In the tables pOIl and |XIII[ we consider a few cases and give a rough estimate 
of the uncertainties due to the errors in fi{q^). The overall error is about ± 5% in the extraction of (^Qs, if we 
consider the errors in fi{q‘^)s are about ± 10%. We also estimate the errors in SCi by considering the actual errors in 
all the form-factors given in ref. |22] , and find that they are even smaller than whatever we have shown in the above 
mentioned tables. Finally, we would like to comment that the estimated errors due to form-factors, and the other 
experimental parameters will have almost equal impact on all the extracted SCiS, which is also small. Therefore, our 
conclusions about the relative sensitivities will not change. 


Cases 


SCi 

<5CC 

SCi- 

±%Err. 

r 

Cl 

0.448 

0.427 

0.472 

5.032 

Re{Cv2 ) 

0.211 

0.201 

0.222 

5.032 

ImlCrY 

0.006 

0.0058 

0.0065 

5.032 

9* 

Cl 

0.042 

0.040 

0.044 

5.032 

Im{Csy 

0.961 

0.916 

1.013 

5.032 

Im{CTY 

0.0048 

0.0046 

0.0051 

5.032 


TABLE XIII: Numerical values of SCiS, and SCC {SCi ) considering +10% (-10%) errors in fi{q‘^) for few cases of R{D*). 

The % error is given by {SCi — SCi^)/SCi. 


IV. SUMMARY 


We use the optimal observable technique to test the sensitivities of the various observables in B ^ to the 

various NP interactions, like new vector, scalar and tensor interactions. Numerically, we find that the observables in 
B —> DrUr are more or less equally sensitive to scalar and tensor interactions, only exception is the r polarisation 
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asymmetry, where 6Re{Cs) < SRe{CT) but they are of same order. Therefore, even if the measured values of the 
observables deviate from their SM expectations, a priori it would be hard to decide what type of new interaction will 
it be. However, the analysis of the distribution of the decay rate allows us to separate the regions of which are 
sensitive to scalar interaction (low and tensor interaction (high g^). 

The overall sensitivity of the observables in H -> D^rUr is more towards tensor interactions, in particular to 
ICtP- Also, we note that 6Re{CT) < SRe{Cv 2 ) but they are of same order, hence, we need better statistics to 
distinguish their effects. The decay B —^ has very poor sensitivity to scalar interaction compared to the tensor 

interaction, the only exceptions being the 1)* polarisation, However, in the absence of tensor interactions, the 

decay B is equally sensitive to real part of both the vector and scalar Wilson coefficients, sensitivity to 

\Cs\‘^ is much less compared to the real parts. However, the analysis of the g^ distributions of the decay rate shows 
distinct regions of g^, which are sensitive to vector, scalar, and tensor interactions respectively. These sensitivities are 
lost in the full g^ integrated observables. Present data on different bins do not have sufficient statistics to conclude 
anything, more precise data could help us to pinpoint the type of NP interaction. Therefore, in an experiment, the 
priority should be given to gaining statistics at those regions of g^. 

We note that both the decay modes are more sensitive to the real part of the coefficients compared to imaginary part. 
Among the various observables, r polarisation asymmetries have better sensitivity to the relevant new coefficients 
{Ci); the uncertainties on the extracted C^s are either less or comparable to that obtained in others. Therefore, future 
data on r polarisation asymmetries could put tighter constraints on the NP parameter space. 
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VI. APPENDIX 

In the following tables, the various /i(g^) for different observables used in our analysis are shown. 


Obs 

A 

R{D) 

tdR 

h 



A 

WhF 

ww 

fs 


80 (l - 

U 



h 

Vr 

Vr 


TABLE XIV: fiS for Rd and r polarisation asymmetry in B ^ Dri^r- 


Obs 

A 

R{D*) 

pR* 

h 


Q* {Rv,+ + Hv- + Hyo) -h 

h 

-20* ((l + |l) {Hlo + 2Hv,+Hv,-) + 

0* ((2 - 'f) {Hlo + 2Hv,+Hv,-) - 

fs 

WHs 

WM 


SG* ((l + + Rt,o)^ 

SG* -(*^T,+ + Rt,- + Rt,o)^ 

h 

^g*jr^HsHv,t 

Vr 

3g*jr^HsHv,t 

Vr 

/e 

-ug* ^ {Ht,oHv,o + Ht.+Hv,+ - Ht,-Hv,-) 

40*^ {Ht.oHv.o + Ht,+Hv,+ - Ht.-Hv.-) 

Vr 

fr 

120. [HTfiHvfi + Ht,+Hv,- - Ht,-Hv,+) 

-40*^ {Ht,oHv,o + Ht,+Hv,- - Ht,-Hv,+) 
Vr 


TABLE XV: fiS for Rd* and r polarisation asymmetry in B ^ D*Th'r- 
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^^Obs 

fi 

aU 


r)R 

/i 

TT f TJS TTS A 

\T* - Hl^) 


/2 

o 

K 


lG*Hl 

h 

V 

1 


sg* (i + ^) Hio 

f4. 


jr*Ji^HsHv,o 

V9 

Vr 

h 

— 

_4jr* {HT,oHv,t + Ht,+Hv,+ + Ht,-Hv,-) 

Vr 

- 12 g*^HT,oHv,o 

Vr 

h 

— 

{HT,oHv,t + Ht,+Hv,- + Ht,-Hv.+) 

Vr 

— 

fr 

- 


- 


TABLE XVI: fiS for r forward-backward asymmetries in B ^ decays, and polarisation asymmetry in 

B D^TVr. 


The expressions for are given by [5] 


g. 

T 


tb 


GllVebf 


B{B Dip) 1927r3m| 


Tb 


G%\Vcbf 2 

ino,..^3 g VAd*(</2) 1 


B{B 1927r^m^ ^ 

tb G%\Vebf 2 
B{B ^ Dip) 128 n^m%‘^ VMq)[^- 

G%\V,bf 2 

100-3^3 g \/Ad.(</2) 1 


Tb 


B{B D'^lu) 1287r^m^ 
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